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Abstract 

It has been shown by A.-L. Barabasi that the priority based scheduling 
rules in single stage queuing systems (QS) generates fat tail behav- 
ior for the tasks waiting time distributions (WTD). Such fat tails are 
due to the waiting times of very low priority tasks which stay un- 
served almost forever as the task priority indices (PI) are "frozen in 
time" (i.e. a task priority is assigned once for all to each incoming 
task). Relaxing the "frozen in time" assumption, this paper studies 
the new dynamic behavior expected when the priority of each incom- 
ing tasks is time-dependent (i.e. "aging mechanisms" are allowed). 
For two class of models, namely 1) a population type model with an 
age structure and 2) a QS with deadlines assigned to the incoming 
tasks which is operated under the "earliest-deadline-first" policy, we 
are able to analytically extract some relevant characteristics of the the 
tasks waiting time distribution. As the aging mechanism ultimately 
assign high priority to any long waiting tasks, fat tails in the WTD 
cannot find their origin in the scheduling rule alone thus showing a 
fundamental difference between the present and the A.-L. Barabasi's 
class of models. 

Keywords : queueing systems - waiting time distributions - fat tails dis- 
tributions - priority indices dynamics - "earliest-deadline first" scheduling - 
tasks with deadlines - age structured population models. 
PACS: 89.75.Da 0.2.50.-r 
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1 Introduction 

In recent contributions A.-L. Barabasi [Barabasi 2005], [Vazquez 2005] and 
[Vazquez et al. 2006] propose a simplified model of the human activity dy- 
namics. These authors view the human activity as a decision based queueing 
system (QS) where tasks to be executed arrive (randomly) and accumulate 
before a server S - here S stands for the processing action of the human 
operator. The time required to process a task (i.e. the service time) is gener- 
ally drawn from a probability distribution. In addition to the usual features 
inherent to any QS, each incoming task is endowed with a priority index 
(PI) which indicates the urgency to process the job. In this conceptual set- 
ting, [Barabasi 2005], [Vazquez 2005] and [Vazquez et al. 2006] study the 
dynamics arising when the service policy is not restricted to the usual first- 
come-first-serve (FCFS) rule but follows scheduling policies based on Pi's. 
Under such priority-based scheduling rules, it is shown that the timing of the 
tasks follows fat tails probability distribution, (i.e the activity of the server 
exhibits bursts separated by long idle periods). This "burst" character has to 
be contrasted with the ubiquitous Poisson behavior which arises when tasks 
are executed according to FCFS or to purely random order scheduling rules. 
In this general context, we shall distinguish between two types of dynamics, 
namely: 

i) Service policies based on fixed (i.e. frozen in time) priority indices. 
This case which is considered in [Barabasi 2005], [Vazquez 2005] and 
[Vazquez et al. 2006] assumes that the value a of the PI is fixed once for 
all. Accordingly, very low priority jobs are likely to never be served. To 
circumvent this difficulty [Barabasi 2005] , [Vazquez 2005] and [Vazquez 
et al. 2006] introduce an ad-hoc probability factor < p < 1 in terms 
of which the limit p — > 1 corresponds to a deterministic scheduling 
strictly based on the Pi's while in the other limit p — > the purely 
random scheduling is in use. In this setting, the waiting time distri- 
bution (WTD) of the tasks before service is shown to asymptotically 
exhibit a fat tail behavior. The main point of the Barabasi's contribu- 
tion is to show that Pi-based scheduling rules can alone generate fat 
tails in the WTD of unprocessed jobs. 

ii) Service policies based on time- dependent priority indices. Here 
the priority index is time-dependent. This typically models situations 
where the urgency to process a task increases with time and a(t) will 
hence be represented by increasing time functions. Clearly, scheduling 
rules based on such time-dependent PI do offer new specific dynamical 
features. They are directly relevant in several contexts such as: 

1) Flexible manufacturing systems with limited ressource. Here a 
single server is conceived to process different types of jobs but only 
a single type can be produced at a given time t (i.e. this is the 
limited ressource constraint). Accordingly, the basic problem is to 
dynamically schedule the production to optimally match random 
demand arrivals for each types of items. The dynamic scheduling 
can be optimally achieved by using time-dependent priority in- 
dices (the Gittins' indices) which specify in real time, which type 
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of production to engage [Dusonchet et al. 2003]. Problems of this 
type belong to the wider class referred as the Multi-Armed Bandit 
Problems in operations research. 

2) Tasks with deadlines. This situation, can be idealized by a 
queueing system where each incoming item has a deadline before 
which it definitely must be processed, [Lehozcky 1996], [Lehozcky 
1997], [Doytchinov ct al. 2001]. In this case, to be later discussed 
in the present paper, we can explicitly derive the lead-time profile 
of the waiting jobs obtained under several scheduling rules, in- 
cluding the (optimal) time-dependent priority rule known as the 
earliest- deadline- first policy. 

3) Waiting time- dependent feedback queueing systems. In queue- 
ing networks, priority indices based on the waiting times can be 
used to schedule the routing through the network. For networks 
with loops, such scheduling policies are able to generate generi- 
cally stable oscillations of the populations contained in the waiting 
room of the queues, [Filliger et al. 2005]. 

In the context of QS, the waiting time probability distribution (WTD), (i.e. 
the time the tasks spend in the queue before being processed) is a central 
quantity to characterize the dynamics. It strongly depends on the arrival and 
service stochastic processes - in particular to the distributions of the inter- 
arrival and service time intervals. The first moments of these distributions, 
enable to define the trafic load p :— — > 0, (i.e. the ratio between the mean 
service time ^ and the mean arrival time j) and clearly the stability of 
elementary QS is ensured when < p < 1. Focusing on the WTD, [Barabasi 
2005], [Vazquez 2005] and [Vazquez et al. 2006] emphasized that heavy tails 
in the WTD can have several origins, three of which are listed below: 

i) the heavy traffic load of the server which induces large "bursty" 
fluctuations in both the WTD and the busy period (BP) of the QS. 
For QS with feedback control driving the dynamics to heavy traffic 
loads, this allows to generate self-organized critical (SOC) dynamics, 
[Blanchard et al. 2004] and the resulting fat tails distribution exhibit 
a decay following a —3/2 exponent. 

ii) the presence of fat tails in the service time distribution produce 
fat tails of the WTD a property which is here independent of the 
scheduling rule [Boxma et al. 2004]. For completeness, we give her a 
short review of these recent results in Appendix A. 

iii) priority index scheduling rules as discussed in [Barabasi 2005], 
[Vazquez 2005] and [Vazquez et al. 2006]. 

Our present paper focuses on the case iii) but contrary to the discussion 
carried in [Barabasi 2005] , [Vazquez 2005] and [Vazquez et al. 2006] , we shall 
here consider the dynamics in presence of age-dependent priority indices. As 
it could have been expected, these aging mechanisms generate new behaviors 
that will be explicitly discussed for two classes of models. 
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2 Scheduling based on time dependent 
priority indices 



The most naive approach to discuss the dynamics of QS with scheduling 
based on time- dependent priority indices is to think of a population model in 
which the members suffer aging mechanisms which ultimately will kill them. 
Naively, we may consider the population of a city in which members are 
either born in the city or immigrate into it at a certain age and finally die 
in the city. Assuming that the death probability depends on each individual 
age, the study of the age structure of the population exhibits some of the 
salient features of our original QS. This is the class of models to be discussed 
in section l2~Tl Later in section l2~2l we shall return to the original model of L. 
Barabasi and consider a simple QS where each task waiting to be processed 
carries a deadline (playing the role of a PI) and as time flows the these 
deadlines steadily reduce - this implies a (time dependence of the PI). At a 
given time, the scheduling of the tasks follows the " earliest-deadline- first" 
(EDF) policy and given a queue length configuration, we shall discuss the 
lead-time (lead-time = deadline - current time) profile of the tasks waiting 
to be served. 

2.1 Tasks population dynamics with time de- 
pendent priority indices 

Consider a population of tasks waiting to be processed by S with the follow- 
ing characteristics: 

i) an inflow of new tasks steadily enters into the queueing system. Each 
tasks is endowed with a priority index (PI) a which indicates its degree 
of urgency to be processed. In general, the tasks are heterogenous as 
the PI are different. In the time interval [t,t + At], the number of 
incoming jobs exhibiting an initial PI in the interval [a, a + Aa] is 
characterized by G(a, t)AtAa. 

ii) Contrary to the situations discussed in [Barabasi 2005], an "ag- 
ing" process directly affects the urgency to process a given task. In 
other words the priority index a is not frozen in time but a = a(t) 
monotonously increases with time t. For an infinitesimal time increase 
At, in the simplest case we shall have a(t + At) = a(t) + At. Here we 
slightly generalize this and allow inhomogeneous aging rates written as 
p(a) > meaning that a(t + At) = a(t) + p(a)At. 

iii) the scheduling policy depends on the PI of the tasks in the queue 
and we will focus on the natural policy "process the highest PI first". 

iv) at time t, a scalar field M(a, t) counts the number of waiting tasks 
with priority index a. Hence M(a,t)Aa is the number with PI p(a) £ 
[a + Aa] . Accordingly, the total workload facing the human server <S 
will be given, at time t by: 




(1) 
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v) in the time interval [t,t + At], the server S processes tasks with 
an a-dependent rate /j(a)At. Typically fi(a) could be a monotonously 
increasing function of a. As the service rate /j,(a) explicitly depend on 
the PI a, it therefore plays an effective role of service discipline. 

The previous elementary hypotheses imply an evolution in the form: 



M (a + p(a)At, t + At)Aa » M(a, t)Aa - fi(a)M(a, t)AaAt + G(a, t)AtAa 

Dividing by AaAt, we end, in the limits Aa — > and At — > 0, with the 
scalar linear field equation: 

d d 

—M(a,t) +p(a)—M{a,t) + fi(a)M(a,t) = G(a,t). (2) 
ot aa 

It is worth to remark that the dynamics given by Eq.(J2J is closely related 

to the famous McKendrick's age structured population dynamics, [Brauer ct 

al. 2001]. 

Assuming stationarity for the incoming flow of tasks (i.e. G(a,t) — G s (a)), 
the linearity of Eq.(J2J) enables to explicitly write its stationary solution as: 



M(a) = 7T (a) 



C 



where 



7r(z) = exp 



G s (z) 



dz 



p(y) 



dy 



(3) 



(4) 



with an integration constant C < oo remaining yet to be determined. Assume 
that the PI attached to the incoming jobs do not exceed a limiting value T, 
namely: 



G(a,t) = l(a < T)G{a,t) 



G s (a)=I(a<T)G s (a), 



(5) 



where I(a < T) is the indicator function. In other words Eq. 101 indicates 
that the new coming jobs do not exhibit arbitrarily high Pi's. 
This enables to define: 



*(T) := / 
Jo 

and Eq.© reads as: 



G s (z) 
p(z)ir(z) 



dz = 



G s {z) 
o p{z)k{z) 



dz 



.,, , J 7r(o) c+ r;-^l 

M(a) = I l Jo p( z )^y z ) 

\ n(a)[C + yj(T)] 



dz 



if 
if 



a<T 
a>T. 



(6) 



(7) 



The asymptotic behavior of M(a) for a — > oo is entirely due to 7r(a), (the 
square bracket terms are bounded by constants) and therefore Eqs. (@J and 
J2} imply: 



f e^r aq when ^4 a fca9_1 wit h q > 0, 

(8) 

In view of Eq. © , the following alternatives occur: 
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a) for q < in Eq.©, the integral J °° M{z) dz does not exist. In this 
case an ever growing population of tasks accumulates in front of the 
server and the queueing process is exploding. 

b) for q > 0, a stationary regime exists and in this case the constant 
C < oo in Eq.JJJ) can be determined by solving: 



which expresses a global balance between the stationary incoming and 
out going flows of tasks. 

c) for q = which implies that oc | , Eq. © produces an exponent- 
Is; fat tail distribution for M{a) counting the number of waiting tasks 
with PI a in the system. For T < oo and a — > oo, the fat tail of M(a) 
is populated by long waiting tasks i.e. those having spent more than 
a — T waiting inside the system before being served. In the limiting 
case, for which /i(a) = /i = const and p(a) = a (i.e. aging directly 
proportional to time) which leads to q — in Eq.©, the density M(a) 
coincides directly with the WTD for a — > oo. 

This population model shares several features with the Barabasi's model, 
namely: 

a) when a stationary regime exists, the function G s (a) which here plays 
the role of the initial PI distribution in [Barabasi 2005] , [Vazquez 2005] 
and [Vazquez et al. 2006], does not affect the tail behavior given by 
Eq.©. 

b) the scheduling rule here is implicitly governed by the service rate 
/i(a) which itself depend on time as the PI a = a(t) are time-dependent. 
Note that /i(a) directly influences the asymptotic behavior of Eq.©. 
In particular for case c), the tail exponent explicitly depends on n(a). 

Besides the similarities, we now also point out the important differences 
between the present population model and the model discussed in [Barabasi 
2005], [Vazquez 2005] and [Vazquez et al. 2006]: 

a) the service is not restricted to a single task at a given time (i.e. the 
service ressource is not limited). Indeed describes an average flow 
of service and hence several tasks can be processed simultaneously - (in 
the city population model the service corresponds to death and several 
individual may die simultaneously). 

b) while the fat tail in [Barabasi 2005], [Vazquez 2005] and [Vazquez 
et al. 2006] is entirely due to the scheduling rule and therefore occurs 
even for QS far from traffic saturation, this is not so in the population 
model. Indeed in this last case, fat tails are due to heavy traffic loads 
occurring when the flow of incoming tasks nearly saturates the server, 
(this is implied by q — in Eq.©) - for lower loads occurring when 
q > the fat tail in Eq.© disappears. 




(9) 
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2.2 Stochastic dynamics. Real-time queueing 
dynamics 

In this section we will use the results of the real-time queueing theory 
(RTQS), pioneered in [Lehozcky 1996], to explore situations where the in- 
coming jobs have a deadline - this problem is already suggested in [Barabasi 
2005]. Based on [Lehozcky 1996], [Lehozcky 1997], [Baldwin et al. 2000] 
and [Doytchinov et al. 2001], first recall the basic hypotheses and the rel- 
evant results of RTQS's. Consider a general single server QS with arrival 
and service being described by independent renewal processes with average 
j respectively — and finite variances for both renewal processes. Each 
incoming task arrives with a random hard time relative deadline T> drawn 
from a PDF G(x) with a density g(x): 

Prob{0 <V<x} = GO), 

with average (D): 

poo roc 

(V)= (1-G(x))dx = xg{x)dx. 
Jo Jo 

At a given time t, we define the lead time C to be given by: 



C = V-t, (10) 

Assume now that the lead time C plays the role of a priority index and the 
service is delivered by using the earliest-deadline-first (EDF) rule with pre- 
emption (i.e. the server always processes the job with the shortest lead time 
C). Preemption implies that whenever an incoming job exhibits a shorter C 
than the one currently in service, this incoming job is processed before, (i.e. 
preempts), the currently engaged task which service is postponed. The EDF 
rule directly corresponds to the deterministic policy (i.e. p = 0, 7 = in the 
original Barabasi's contribution [Barabasi 2005]. 

At a given time, one can define a probability distribution corresponding to 
the lead time profile (LTP), Prob{— 00 < C < x} := F(x), of the jobs waiting 
in the QS. The LTP specifies the repartition of tasks having a given C at 
time t. Knowing the queueing population Q at a given time, it is shown 
in [Doytchinov et al. 2001] that for heavy traffic regimes, the LTP can, in 
a first order approximation scheme, expressed by a simple analytical form. 
Specifically, following [Doytchinov et al. 2001], define a frontier J-(Q) > 
to be the unique solution of the equation: 



9. 

A 



/ (1 - G(x)) dx, and G{x) = for x < 0). (11) 

I HQ) 



In [Doytchinov et al. 2001], it is shown that two alternative regimes can 
occur: 

a) Jobs served before deadline. Solving Ea. (|llf) . a positive value J~(Q) > 
is obtained. In this case, the LTP cumulative distribution F(x) takes 
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the form, see Figure ^ 



F(x) 



when x < T(Q), 

A (J~ [1 - G( V )j drf) when T(Q) < x. 



(12) 



b) Jobs served after deadline. Solving Ea. pi(l . a negative value J-(Q) < 
results. In this case, the LTP cumulative distribution F(x) takes the 
form, see Figure 



r o 



F(x) = < 



x-(V) + % 



where q = 



CD}' 



when x < T(Q) = (V) - f < 0, 
when F(Q) < x < 0, 



(13) 




Figure 1: Density of the lead time profile f(x) = dF }^ when T{Q) > 0. 




FIX) *ffe) 



Figure 2: Density of the lead time profile fix) = ^j^- when J-(Q) < 0. 



Remark. The alternative regimes given by Eas. (|12fl and i|13|) can be heuris- 
tically understood by invoking the Little law which connects the average 
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queue length (Q) with the average waiting time (W), [Cohen 73]: 



(Q) = X(W), (14) 

which is independent of the scheduling policy. In view of Eqs.^TJ and i|14() . 
one obviously suspects that the LTP strongly depends on the sign of the 
difference (D) - = (V) - (W). Intuitively, when (W) exceeds (D), it is 
expected, in the average, that processed jobs will be delivered too late and 
conversely. While the above heuristic arguments is strictly valid only for the 
averages, [Doytchinov et al. 2001], were able to show that in heavy traffic 
regimes, it also holds also for the LTP given in Eos. (|12(l and 



Assuming that the arriving tasks have positive deadlines, the LTP as given 
by Eqs.JEI) and {T3J| imply: 

a) If the left-hand support of the LTP is negative, then a job entering 
into service is already late, (case of Ea. (|13|l % l see Figure [21 

b) If the left-hand support of the LTP is positive then a jobs enters 
into service with a positive lead time, (case of Eq. ((T2T ) see Figure 
Accordingly, it is likely that the tasks will be completed before the 
deadline expired. 

c) The critical value Q* = for which J-(Q) = 0, corresponds to a 
queue length for which customers are likely to become late. Choosing 
Q exactly to Q* , we cannot expect lateness to disappear completely 
but for Q < Q* lateness will be strongly reduced a behavior clearly 
confirmed by simulation experiments [Doytchinov et al. 2001] and 
[Baldwin et al. 2000]. 

d) For deadline distributions G(x) with fat tails, it follows immediately 
from Eas. H12|> and ljl3|l that the LTP does possess a fat tail. 



2.2.1 "First come first served" (FCFS) scheduling poli- 



cies 



For the choice g(x) = 5(x), (i.e. zero deadline), the EDF scheduling policy 
directly coincides with the FCFS rule. Indeed, in this case T(Q) = and the 

LTP density is given by Eq. Ijl2(l is a uniform probability density U 



Q 



_w 



being its support). This expresses the fact that in the heavy traf- 
fic regime p — A//i ~ 1, the waiting time behaves as Q x (i) « Q x (j) 
leading to a LPT linearly growing with Q. For general G(x) 1 the LTP asso- 
ciated with a FCFS scheduling rule will be given by the convolution of the 



_Q o 



Indeed, 



deadline distribution G(x) with the uniform distribution U 

adding the task deadlines with the time spent in t5he queue, we recover the 
tasks lead-time. Accordingly, in the heavy traffic regime and for a given 
queue length Q, one explicitly knows the LTP's for both the EDF and the 
FCFS scheduling policies thus enabling to explicitly appreciate their respec- 
tive characteristics. In particular, using Eas. (|12(l and (|1 H|l . one can conclude 
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that for a given queue length Q the FCFS scheduling rule the associated LTP 



F{x) being the convolution of G(x) with the U 



.9. o 



, it takes the form: 



F(x) 



A f x 

Q J- 



f ) * 



when x < — 



when - ^ < a; < 0, 



. «+^{r[G(e+§)-G(o]^} 

where the constant k reads as: 



when x > 0, 



(15) 



A 

Q ./_(«) 



G(C + f) 



Ea. fl5|) allows to emphasize the following features: 

i) When the left-hand support of the deadline distribution G(x) is 
larger than y , the left boundary of the support of F(x) is larger than 
and therefore the jobs experience no delay when entering into service. 

ii) If the left-hand support of G(x) is smaller than then it may 
happen that the LTP exhibits a negative left-hand support under the 
FCFS policy and a positive left-hand support under the EDF schedul- 
ing rule. Hence in this last situation, the FCFS policy would deliver 
tasks with lateness while the EDF tasks will be processed in due time. 
This explicitly confirms intuition that EDF is indeed an efficient pol- 
icy. It has been shown that the EDF scheduling rule is optimal for 
minimizing the number of jobs processed after the deadline [Panwar et 
al. 881. 



iii) If G{x) exhibits a fat tail for x — » oo so has the LTP and this 
whatever the scheduling rule in use. This can e directly verified from 
Eq.JEIJ by studying the LTP density f(x) = for x 

have: 



oo, we 



/(*) 



X 

Q 



G(x + t) - G(x) 
A 



for 



x — > oo, 



which when G{x) ~ 1 — x q and for y < const takes the form 



fix) 



"(9+1) 



for 



(16) 



Hence, the LTP inherits the fat tail property of G(x) and this even 
when using the optimal EDF scheduling rule - a fully explicit illus- 
tration involving the Pareto probability distribution is given in the 
Appendix B. 

The results obtained for the LTP, enable to get the asymptotic properties of 
the waiting time distribution WTD. Indeed, under the EDF policy, the more 
urgent jobs are served first and therefore the waiting time before service of 
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the queueing jobs will be larger or equal to the lead-time. Accordingly, if the 
LTP exhibits a fat tail distribution so will the WTD. Hence, while the EDF 
policy decreases, compared with the FCFS rule, the number of jobs served 
after their deadline, it cannot get rid of the fat tails generated by the deadline 
probability distribution G(x). Let us emphasize here, the fat tails of the LTP 
(and hence of the WTD) are here entirely due to G(x) and his asymptotic 
behavior of the LTP is shared by both the EDF and FCFS policies. This 
is fundamentally different from the frozen in time PI models discussed in 
[Barabasi 2005], [Vazquez 2005] and [Vazquez et al. 2006] where the fat tail 
behavior does not depend on G(x). This can be heuristically understood as, 
in [Barabasi 2005], [Vazquez 2005] and [Vazquez et al. 2006], the fat tail is 
mainly due to the low priority jobs which, as no aging mechanism occurs, 
are likely to never be served. Note that in [Barabasi 2005], [Vazquez 2005] 



and [Vazquez et al. 2006], stable queueing models, (i.e. those for which the 
traffic p < 1), fat tails of the WTD disappear under a FCFS scheduling rule. 
Indeed without priority scheduling, the WTD always follows an exponential 
asymptotic decaying behavior. In presence of time dependent PI, all tasks 
do finally acquire a high priority and this aging mechanism precludes the 
formation of a fat tail solely due to the scheduling rule. Accordingly, in 
presence of aging PI, the generation of WTD with fat tails will be due to 



There are several possibilities to analytically discuss the scheduling of tasks 
in QS with time dependent priority indices and to infer on the existence of 
fat tails for the asymptotic behavior of the resulting WTD. In this note, we 
propose two distinct models where an explicit analysis can be developed. 
Our first model is directly inspired by the study of age classes in population 
dynamics where the mortality rate increases with the age of the individuals. 
In this context, identifying the service of the QS with the death of an indi- 
vidual, this dynamics is closely related to the scheduling based on PI, the 
indices here being the age of the individuals and the immigration with dif- 
ferent ages plays the role of incoming tasks with different priorities. For this 
class of dynamics, it is straightforward to show that fat tails of the WTD can 
develop on the onset of stability of the population model. As in the original 
Barabasi model, the tail behavior of the WTD does not depend on the de- 
tail nature of the PI but only on the scheduling rule - (corresponding in the 
population model to the mortality rate). Our second modeling frame which 
is closer to the Barabasi original idea, we consider a classic QS in which the 
scheduling rule is based on the deadlines attached to each incoming tasks. 
As time flows, the deadlines reduce and hence the waiting tasks acquire a 
higher priority to be processed. In the heavy traffic limit, i.e. for regimes 
where the law of large numbers dominate, it is possible to analytically derive 
the lead-time profile (lead-time = deadline minus the time elapsed in queue- 
ing) of the waiting tasks and from this to get information on the asymptotic 
behavior of the associated WTD. In this case and contrary to the conclu- 
sions exposed in [Barabasi 2005], [Vazquez 2005] and [Vazquez et al. 2006], 



G(x). 



3 Conclusion and summary 



n 



the scheduling rule solely cannot generate fat tails in the WTD. Fat tail in 
[Barabasi 2005], [Vazquez 2005] and [Vazquez et al. 2006] are due to low 
priority jobs which are likely to never be served. This possibility disappears 
when time-dependent PI are considered as, due to aging, initially low pri- 
ority tasks do acquire, with time, high priorities and hence will not stay 
unprocessed forever. This precludes the formation of fat tails in the WTD. 
We finally observe that, in this second class of models, the only possibility 
to generate fat tails is to feed the QS with tasks deadlines drawn from a fat 
tail distribution. 



Acknowledgements. M.O.H. thanks numerous fruitful discussions with 
Olivier Gallay and Roger Filliger. 
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Appendix A - Waiting time distribu- 
tions for QS with fat tail service times 



Let us reproduce here a result recently derived in [Boxma et al. 2004]: 

Theorem 1 . Assume that the (random) service time in a M/G/l QS is 
drawn from a PDF with a regularly varying tail at infinity with index v G 
(— 1, —2), (regularly varying with index v G (—1, —2) =>■ fat tail with index 
v G ( — 1, —2)). For this range of asymptotic behaviors of the PDF, the first 
moment (3 of the service exists. Assume further that the service is delivered 
according to a random order discipline. Then the waiting time distribution 
W'ros exhibits a fat tail with index (1 — v) € ( — 1, 0) and more precisely, we 
can write: 

Prob(^ ROS > x) a C-^- r x l ~ v C{x), (17) 

1 - p (3 T(2 - v) 

where p < 1 is the traffic intensity, /3 the average service time, C(x) a slowly 
varying function and 

h{p,v) ■= I f(u,p,u)du, 

with: 



The fat tail behavior given in Ea. H17|l is therefore entirely inherited from 
the fat tail behavior of the service and is not affected by any reduction of 
the trafic intensity p. Note also that change of the scheduling rule cannot 
get rid of this fat tail behavior. This point can be explicitly observed in 
[Cohen 1973], [Pakes 1975] who show that for the previous M/G/l QS with 
a random order service (ROS) service discipline, one obtains: 

Prob(WRos > x) cc h(p, ^)Prob(WFCFS > x), for x — > oo, (18) 

from which we directly observe that the fat tail in the asymptotic behavior 
in not altered by a change of the scheduling rule. 

Note finally that for the M/M/l QS, (i.e. exponential service distributions 
and hence no fat tail), [Flatto 1997] shows that the random order service 
scheduling rule yields: 

Prob(W RO s > x) a C p x~i e -^~ Sxi , for x -> oo, (19) 

with 
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(7(p) = 23 3 2 7r6pi2- — gexp> 



7 = 



(,-» 



- 1 



1 — ^2 



and (5 = 3 



1 — pa 



which has to be compared with the FCFS scheduling discipline, which for 
the same M/M/l QS reads as, [Cohen 1973]: 



Prob(WW s > a:) = ~(1 - p)e-^- p >. 



(20) 



While the detailed behaviors given by Ea. H19(l and (|20|) clearlv differ, they 
however both share, in accord with [Barabasi 2005], an exponential decay. 



Appendix B - Deadline drawn from Pareto 
distribution 



Here, we focus on: 



G(x) 



when ■§ < 1, 



l-(f) 



(21) 



when ■§ > 1, oj > 1, 



which has no moment of order u — 1 or higher. For to > 2, we have (I?) = 



£: 



5 



BX 



Q(uj-2) 



u>-2 i -° A 



when ^ < -^>, 



when ^ > — ~~n . 



(22) 



Using Eas. (|12|) and JT3J, the LTP distribution reads as: 



I? 



— > 

\ ~ u - 2 



' 



x {u) ~ x B-x 



1 - 



BX IB 



(f) D 



when x < T{Q), 
when ^(Q) < x < B, 



when x > B. 



(23) 



Q < B 



X uj-2 



F(x) 



1 - 



BX ( B 



Q(w-2) 



(f) 



when x < J-(Q), 



when x > J-{Q). 



(24) 
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Eqs.(H3J) and JUJl exhibit a fat tail with power uj — 2. Note that, Ea. lf^|) . 
implies that for to > 2 and for y < rjzm , the EDF scheduling policy part 
of the tasks enter into the service before the due date expired. Finally note 
also, that for ui < 2, no moments exists for the deadline distribution and 
hence the theory [Doytchinov et al. 2001] cannot be applied directly. We 
conjecture that for these regimes no scheduling rule will be able to deliver 
tasks in due time. 
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